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Abstract 

We study polynomials that are orthogonal with respect to a varying 
quartic weight exp{—N(x'^/2-\-tx'^/4)) for t < 0, where the orthogonal- 
ity takes place on certain contours in the complex plane. Inspired by 
developments in 2D quantum gravity, Fokas, Its, and Kitaev, showed 
that there exists a critical value for t around which the asymptotics 
of the recurrence coefficients are described in terms of exactly spec- 
ified solutions of the Painleve I equation. In this paper, we present 
an alternative and more direct proof of this result by means of the 
Deift/Zhou steepest descent analysis of the Riemann-Hilbert problem 
associated with the polynomials. Moreover, we extend the analysis 
to non-symmetric combinations of contours. Special features in the 
steepest descent analysis are a modified equililbrium problem and the 
use of ^'-functions for the Painleve I equation in the construction of 
the local parametrix. 
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1 Introduction and statement of results 



1.1 Introduction 

Let V = Vt he the quartic polynomial defined by 

V{x) = Vt{x) = xV2 + txV4, (1.1) 

where t G M is a parameter. In this paper we study orthogonal polynomials 
with respect to the exponential weig ht e-^^*(^) where N eN, in cases where 
t is negative. Since for negative t integrals like 

are divergent, we will have to define what we mean by orthogonal polynomials 
in this case. Before doing so, and in order to motivate what we are going to 
do, we discuss the case t > first. 

For t > 0, the monic orthogonal polynomial 7r„ of degree n G N satisfies 

7r„(x)x'= e-^^*(")dx = 0, (1.2) 

for A; = 0, 1, . . . , n — 1. An important feature of these polynomials is that the 
recurrence coefficients a„ in the three term recurrence nn+i{x) = xiinix) — 
ttnTin-iix) satisfied by these polynomials satisfy 

Tl 

a„ + ta„(a„_i + a„ + a„+i) = — . (1.3) 

The nonlinear difference equation ()1.3p is referred to as the string equation 
or the Freud equation. It is also known as an example of a discrete Painleve 
equation, see j27j and the references cited therein. Of course the polynomials 
7r„ as well as the recurrence coefficients a„ depend on and t. If we want 
to emphasize this dependence we write nn,N,t and an,N(t)- 

An important tool in the study of the asymptotic behavior of Tin,n,t and 
0'n,nit) as n — > cxD is played by the equilibrium measure in the external field 
Vt [211 ■ This is the unique Borel probability measure /i = /if on R that 
minimizes 

Iv^{^l)= //log A ^i{x)d^i{y) + Vt{x)d^i{x) (1.4) 
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among all Borel probability measures on M. The measure /Xi can be calculated 
explicitly. It is supported on the interval [— q, q] where 



c; 



(1,5) 



and it has a density given by 



dx 



t 



for X e [-Ct, Ct] 



(1.6) 



where 




(1,7) 



The formula p.7p may look a bit awkward since < for t > 0, but we 
have written it this way, since we will mainly use ()1.7|) for t < and then 



The limiting behavior of the recurrence coefficients a„^„(t) as n — > cxd is 
directly related to the support [— Cf,Cf] of the equilibrium measure, since we 
have for t > that 



as n ^ oo. The asymptotics ()1.8|) follows from far more general results in fH] 
where it is also shown that 0{l/n) term can be written as a full asymptotic 
expansion in powers of 1/n. 

1.2 Critical behavior of recurrence coefficients 

As said before, we are going to consider the case t < 0. Although in this 
case the orthogonal polynomials are not well-defined by (jl.2p . the Freud 
equation (jl.3p makes perfect sense. If t > —1/12 also the measure (|1.6|) 
with Ct > and dt > given by ()1.5|) and ()1.7|) is well-defined and gives 
a probability measure on M. This measure does not minimize ()1.4j) among 
all Borel probability measure on M (in fact there is no such minimizer), but 
it does minimize (jl.4p among all Borel probability measures on [— Cf,Q], or 
among all Borel probability measures on the larger interval [—dt,dt\. The 
value t = —1/12 is critical, since for t < —1/12 the measure ()1.6p is not 



d^ > 0. 



VI + 12t - 1 
6t 



(1.8) 
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tcr = —1/12, we have = = 8 
Her where 

for a; G [-V8, VS]. (1.9) 

By formal calculations based on the Freud equation (|1.3|) it was shown 
in the physics literature, see e.g. the surveys [1^1 CHI, that the recurrence 
coefficients an,n{t) have very interesting limit behavior as — > cxd for t near 
the critical value tcr = —1/12. Namely, if t depends on n and tends to tcr as 
n — i> cxD in such a way that 

n^^^{t-tcr) = -cix (1.10) 

remains fixed, then 

lim n2/5(a„_„(t) - 2) = C2y{x) (1.11) 

n^oo 

where y{x) is a solution of the Painleve I equation 

y" = 6y^ + x. (1.12) 

In ()1.1()|1 and p.llj) we have that ci and C2 are certain explicit positive 
constants. 

In two very important papers jTHl Fokas, Its, and Kitaev were able 
to prove this result in a mathematically rigorous way. First of all they made 
it clear how the orthogonal polynomials with respect to Q-^^ti^) should be 
defined in case t < 0. The solution is to consider the orthogonality not on 
M but on a contour F in the complex plane chosen so that Re Vt{z) +oo 
as z ^ oo along the contour. For symmetry reasons they chose a contour 
F that consists of the two lines parametrized hy z = r^™!'^ and z = re~*'^/^ 
with — oo < r < oo. Then the orthogonality property for the polynomials 7r„ 
is 

J 7r„(z)2V^^*('M2 = 0, for = 0, 1, . . . , n - 1, (1.13) 
and the integrals are well-defined. Note that the bilinear form 

{p,q) = j piz)qiz)e-^^'^'Mz (1.14) 



well-defined anymore. Note that for t - 
and we find the critical measure fi-i/u 



d^^^(g_ 2)3/2 

dx 2471^ ^ ' 
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is not positive definite, so that (jl.l7|) is an example of non-Hermitian ortliog- 
onality. As a result it is not automatic tliat a unique monic polynomial 7r„ of 
degree n exists that satisfies p.l3|) . However, it follows from the analysis of 
[THl [T^ that in the asymptotic regime given by (jl.lOj) the monic polynomi- 
als Tin+i,n,t, '^n,n,t, '^n-i,n,t exist for u large enough, they satisfy a three-term 
recurrence relation with recurrence coefficient a„,„(t) that satisfies (jl.llj) for 
a uniquely specified solution y{x) of the Painleve I equation (jl.l2j) . 

The analysis in ^Hj is based on a characterization of the orthogonal 
polynomials by a Riemann-Hilbert problem. The Riemann-Hilbert problem 
admits a Lax pair formulation, and the Freud equation for the recurrence 
coefficients follows from the compability equation for the Lax pair. In [T2] the 
Painleve I asymptotic behavior is proved via a WKB analysis of the Riemann- 
Hilbert problem and the associated Lax pair. It is interesting to note that 
the Riemann-Hilbert problem for orthogonal polynomials was stated for the 
first time in (TH]. It should also be noted that [12] is not restricted to the 
quartic potential p.ll) but more general polynomial potentials are considered 
as well. 

Subsequent developments in the asymptotic analysis of Riemann-Hilbert 
problems made us consider this problem again. The main new development 
is the application of the powerful Deift/Zhou steepest descent analysis for 
Riemann-Hilbert problems |T3] to the Riemann-Hilbert problem for orthog- 
onal polynomials by Deift et al. ^21 which had a major impact on the 
asymptotic theory of orthogonal polynomials. 

In this paper we give an alternative proof of the Painleve I asymptotics 
using these new techniques. The proof we thus obtain is a direct proof, in 
the sense that the analysis only deals with the Riemann-Hilbert problem 
and not the Lax pair. We derive asymptotics for the recurrence coefficients 
without considering the Freud equation in an explicit way. For clarity of 
presentation we restrict ourselves to the quartic potential fjl.H) but extension 
to more general potentials is possible. But even in the context of the quartic 
potential our method allows an extension to a more general model which we 
describe now. 
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Figure 1: The contour T with the weights on the different parts of F 



1.3 Statement of results 

The orthogonality relation with parameters a and (3 

Instead of the bilinear form (jl.l4|) we put weights on the different parts 
of the contour F and consider a bilinear form 

{p,q) = Y^a, p(z)g(z)e-^^'(")dz (1.15) 
i=i "^r, 

where Fj is the part of F in the zth quadrant. Not all choices of weights are 
relevant. We impose ai + 014 = 1 and 02 + 03 = 1. Then putting a = ai 
and /3 = 03 we obtain the situation indicated in Figure ^ The bilinear form 
then depends on the two parameters a and [3 and we denote it by (p, q)a,i3 
so that 

(p,g)„^^ = a I p{z)q{z)e~^^'^''^<lz + {I - 13) I p{z)q{z)e-^^'^''^<lz 

+ (3 [ p{z)q{z)e-^^'^'^dz + {l-a) [ p{z)q{z)e-''^'^'Mz. 

(1.16) 

The choice a = (3 corresponds to the situation considered by Fokas, Its, and 
Kitaev 

The orthogonal polynomial 7r„ of degree n is now defined by the relations 
(7r„(^), z^)^^p = 0, for A; = 0, 1, . . . , n - 1. (1.17) 
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The polynomial clearly depends on a and (5. It is again an example of non- 
Hermitian orthogonality, so that uniqueness and existence are not immediate. 
However if three consecutive monic polynomials 7r„+i, 7r„ and 7r„_i exist then 
they are connected by a three-term recurrence relation of the form 

vr„+i(2) = {z- 6„)7r„(2;) - a„7r„_i(z) (1.18) 

where we now have two recurrence coefficients a„ and hn- The 6„ coefficient 
vanishes only if a = /3, since then the bilinear form is even. In that case 
a„ satisfies the Freud equation ()1.3|) . For a ^ [3 the bilinear form is not 
even and we have the general three-term recurrence (|1.18p . Then (|1.3|) is 
not satisfied, but instead the a„ and satisfy a more complicated system of 
nonlinear difference equations. 

We use an,N{f) and bn,Nit) to indicate the dependence on N and t. The 
recurrence coefficients also depend on a and j3 but we will not indicate that 
in the notation. 

The regular case -1/12 <t<0 

We refer to the case where t is fixed with —1/12 < t < as the regular 
case. In the regular case we will find that the asymptotic formula for a„^„(t) 
is a straightforward continuation of the one valid for t > 0, see (II. 8|) . In 
addition, the recurrence coefficient &„,n(^) is exponentially small. This is our 
first result. 

Theorem 1.1. Let — l/12<t<0 be fixed. Then there exists an uq = no{t) 
such that the recurrence coefficients an,n{t) and bn,n{t) exist for all n > no, 
and 

anAt) = ^^^^^^ + 0{l/n), n^oo, (1.19) 
or 

and there exists a constant Ct > such that 

bnAt) = C(exp(-Ctn)), n ^ oo. (1.20) 

The recurrence coefficient a„_„(t) has a full asymptotic expansion in powers of 
and all terms in the asymptotic expansion are independent of the choice 
of a and /?. 

The proof of Theorem II. II is given in Section |21 
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Figure 2: Jump matrices for the ^'-function associated to the solution ya of 
the Painleve I equation 

The critical case 

The main purpose of the present paper is to deal with the critical case where 
t depends on n such that t— >— l/12as?2— >ooata critical rate. It is only in 
the critical case that the dependence on the parameters a and j3 plays a role. 
In our main result (Theorem II .21 below) we will give asymptotic formulas for 
o.n,n{i) and bn,nit) as n — >• cxD in terms of special solutions Ua and 1/13 of the 
Painleve I equation ()1.12|) . We describe these special solutions first. 

Solutions Ha of the Painleve I equation 

To explain how the solution i/a depends on a we recall the Riemann-Hilbert 
problem associated with Painleve I. Let \E'(C;x) be a 2 x 2 matrix valued 
function that is analytic in the complex ^-plane except for jumps on the 
contours a.Tg( = ±27r/5, arg^ = ±47r/5, and a.Yg( = n. The contours are 
oriented as in Figure |21 The orientation induces a + and — side on each part 
of the contour, where the + side (— side) is on the left (right) if one traverses 
the contour according to the orientation. For ( on the contour the limiting 
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values $±(C; from both sides exist and are connected by 

vI/+(C;x) = vl/_(C;a;)A„ 

where the jump matrices Aa are given in Figure |21 Furthermore, \I/ has 
asymptotics 



for ^ — »• cxo, where o"3 = ( q -i ) ^-^d 9 is defined by 

^(C, a:) = ^C'/' + (1.22) 
5 

Then it is known that 

l/„(x) = 2i(vl/2(x))i2 (1.23) 

is a solution of the Painleve 1 equation, and it is this solution tja that will 
appear in our main result below. 

As is the case for any solution of the Painleve I equation we have that t/a 
is a meromorphic function with an infinite number of poles in the complex 
plane. The Riemann-Hilbert problem for \E' has a solution if and only if x is 
not a pole of i/a- 

The above Riemann-Hilbert problem does not generate all solutions of 
the Painleve I equation, but only those solutions that satisfy 



y(x) = ^/-x/6{l + o{l)) asx^-oo. (1.24) 
All ya satisfy p.24|l and they have a common asymptotic series 



1 + ^^.r_^^-5fc/2 



ak[-x) 

k=l 



as X — ^ — oo (1-25) 



for certain coefficients a^. 

For the special value a = 1 we have that the asymptotic series p.25|) is 
valid as |x| — oo with argx G [37r/5,7r]. The other solutions differ from yi 
by exponentially small terms only. It follows from results of Kapaev [221 ^^^^ 

as X — »• — oo. The behavior (jl.26|) characterizes i/a- 
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Main result 

Our main result is the following theorem. 

Theorem 1.2. Let G C, and let t vary with n such that 

n^l\t + 1/12) = -cix, ci = 2-9/53"'^/^ (1.27) 

remains fixed, where x is not a pole of Ua and yp. Then, for large enough n, 
the recurrence coefficients an^nif) and hn,n{t) associated with the orthogonal 
polynomials with respect to the bilinear form fll.l6|) exist and they satisfy 

a„,„(t) = 2 - C2 (i/,(x) + y^{x)) n^^l'' + 0{n-^"'), = 2^/''?>'"\ (1.28) 

and 

bnAt) = C3 Mx) - y^{x)) n-'/' + 0{n-'/'), c, = 2'/''3'/', (1.29) 

as n ^ oo. The expansions (jl.28|) and (jl.29|) hold uniformly for x in compact 
subsets of M not containing any of the poles of ya and yjs, and the O terms 
can be expanded into a full asymptotic expansion in powers of n~^/^ . 

Remark 1.3. It is quite remarkable that we do not have a term of order 
^-1/5 in (jl.28j) and (|1.29|) . The terms in these expansions ultimately come 
from the terms in the asymptotic expansion (|1.2H) of \E'(C; a;), as we will show. 
Terms of order n~^^^ could have appeared because of the term \l'i(a;)/C^/^ 
in p.21|) . However, it turns out that the entries of '^i{x) cancel out in the 
calculations. These entries also do not contribute to the n~'^/^ term. 

Remark 1.4. Not much is known about the precise location of the poles 
of ya on the real line, however see |Hl 1^- Because of (jl.24|) we know that 
there can be only a finite number of poles on the negative real axis. Joshi 
and Kitaev Prop. 3] showed that every real valued solution of ()1.12p has 
at least one pole on the positive real axis, so this applies in particular to our 
solutions ya with Re a = 1/2. 

Remark 1.5. Even though we will concentrate on the quartic potential, it 
will be obvious that our method generalizes to higher degree potentials. The 
quartic potential serves as a generic example for all polynomial potentials 
for which the density of an associated critical equilibrium measure vanishes 
with an exponent 3/2 at the endpoints of the support as in (jl.9p . 
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Note that this type of vanishing is not possible in the case of usual or- 
thogonahty with respect to varying exponential weights on the real line, since 
then it is only possible that the density of the equilibrium measure vanishes 
at an endpoint with an exponent {Ak + l)/2 with k G Nq. See |1U1 [TT| for 
more details. The generic case k = leads to Airy functions. The first criti- 
cal case = 1 is described in terms of the second member of the Painleve I 
hierarchy j7j. 

Remark 1.6. The general symmetric quartic potential V{x) = gx^ / 4: 
exhibits another type of critical behavior for (7 < and t > 0. Here the zeros 
of the orthogonal polynomial 7r„ with respect to e~"^^^^ may accumulate 
either on one or on two intervals, depending on the values of g and t, and the 
same holds true for the eigenvalues of the unitary random matrix ensemble 
(l/Z„)e-"'^''^(^)rfM on n X n Hermitian matrices M. Bleher and Its |2l El 
showed that the transition is described by the Hastings-Mcleod solution of 
the Painleve II equation. This result was generalized in |B] to more general 
potentials. 

Overview of the rest of the paper 

The proofs of Theorems 11.11 and 11.21 are based on a steepest descent analysis 
of the RH problem that characterizes the orthogonal polynomial 7r„. This 
RH problem is due to Fokas, Its, and Kitaev JJ^. We discuss it in detail in 
Section |21 

In Section El we prove Theorem 11.11 After the appropriate deformation of 
contours, the steepest descent analysis is similar to the analysis in Deift et 
al. IJJJ and therefore we will not give all the details of the proof here. 

The major part of the paper is devoted to the proof of Theorem II . 21 which 
is given in Section 01 Here the critical measure (jl.9|) corresponding to the 
critical value t = —1/12 comes into play. For t 7^ —1/12, we will not use 
the equilibrium measure fit from ()1.6|) in the steepest descent analysis, but 
instead a modified equilibrium measure ut (in fact a signed measure). Similar 
modified equilibrium measures were used before in jH EHH] for double scaling 
limits arising in critical random matrix models. The modified equilibrium 
problem is discussed in Section 4.2. It leads to a modified (yf-function [T^ ITT] 
which is used in the first transformation of the RH problem. 

The steepest descent analysis then proceeds as in [TT]. We open lenses 
around the critical interval [— \/8, \/8] and then construct a parametrix for the 
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Figure 3: Deformation of the contour 



resulting RH problem. The \l/-functions associated with Painleve I appear 
in the construction of the local parametrices around the endpoints ±v^. 
This construction is analogous to the construction in jSJ E], where a local 
parametrix was built out of the \E'-functions associated with Painleve II. The 
final transformation of the RH problem is given in Section 4.8. It leads to a 
RH problem which for large n can be solved explicitly in a series involving 
powers of n~^/^. We calculate the terms up to order n"^/^ explicitly. We 
also show how the recurrence coefficients can be expressed in terms of the 
final RH problem and then the final computations in Section 4.11 lead to the 
proof of Theorem 11.21 

2 The Riemann-Hilbert problem 

By analyticity, a deformation of the contour F does not change the polyno- 
mial 7r„ as defined in (|1.17|) . provided that each part F^ extends to infinity 
in the sector 

Sj = {zeC\ (4j - 3)7r/8 < arg(^) < (4j - l)7r/8}. (2.1) 

Figure 01 shows a particular deformation of F, which will turn out to be 
convenient for our analysis. The contours are deformed so that they contain 
the part [— c?t, on the real axis where dt is defined by fll.7|) . At the points 
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±dt the contours separate and extend to infinity in the sectors Sj. The 
precise setting of this deformation will be given later. Changing our previous 
notation, we will now use Fi, . . . , r4 to denote the contours as in Figure El 
We also put Fq = [—dt, dt]. The contours are oriented as indicated in Figure 

lU 

The parts Fi, . . . ,F4 carry the weights as before, and Fq has the weight 
1. If then for z G F we define 



r 1, z e Fo, 

a, z E Fi, 
1-/3, 2GF2, 
ze F3, 
1 — a, z E F4, 



then for polynomials p and g, the bilinear form p.l5|) is equal to 



o,/3 



p{z)q{z)a{z)e~^^'^''^Az. 



(2.2) 



Consider the following Riemann-Hilbert (RH) problem for a 2 x 2 matrix- 
valued function F : C \ F C^^^^ 



Y{z) is analytic in C \ F 



YAz) 



Y^{z) 



1 a(2)e-^^*(^) 







Y{z) = {I + 0{l/z)) 



1 

n Q 



2 e F 



(2.3) 







00, 



By a standard argument one can show that if a solution to the RH problem 
(j2.3|) exists, it is unique, see [121 ■ The existence of Y depends on the 
existence of polynomials that are orthogonal with respect to the bilinear 
form ()2.2|) . To be precise, we have 

Proposition 2.1. The RH problem (j2.3p has a solution if and only if there 
exist polynomials p and q such that 



(a) p is monic of degree n and {p{z), z'')a,i3 = for k = 0,1, 



.n- 1, 



(b) q has degree < n — 1, {q{z), z^)a,i3 = for k = 0,1, ... ,n — 2, and 
{qiz),z^-')^,f3 = -2m. 
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In that case the solution of the RH problem is given by 




1 f pis)a{s)e-^^ \ 



1 f g(s)a(s)e-^^*W 



zeC\T. 



(2.4) 



■ S 



Tri /r s — z 



Proof. The proof is standard, see JH] or [H]. 



□ 



Note that the non-Hermitian orthogonahty is exactly the right concept 
for the formulation of a RH problem. This has been exploited before in for 
example [H 121 123- 

It follows from Proposition 12 . II that if Y exists, then p = Yn = 7r„ and so 
the monic orthogonal polynomial 7r„ exist. The polynomial q = Y21 satisfies 
the orthogonality conditions for 7r„_i. If q has degree n — 1, then the monic 
orthogonal polynomial 7r„_i exist, and there is a constant Hn-i 7^ such that 
q = Y21 = Kn-iT^n-i- It might happen that q has degree < n — 1, and then 
7r„_i does not exist. 

We use 

y(n+i)^ to denote the solutions of the RH problem 

()2.3|) for the values n + 1, n, and n — 1, respectively (with fixed A^). If all 
three RH problems have a solution then Hn+i, tt^ and 7r„_i exist. We show 
that they are related by a three term recurrence relation. 

Proposition 2.2. Suppose that the monic orthogonal polynomials iTn+i, T^n, 
7r„_i exist. Then they satisfy 



for certain coefficients an and 6„. 

Proof. From Proposition 1)2.11) it follows that the RH problems for F^"^ and 
are solvable. From the discussion above we find that the first column 



Since the jump matrices in the RH problem ()2.3|) have determinant one, 
and detF*^"'^(z) — > 1 as z — 00, it is easy to check that detF'^"^(z) = 1 for 
all 2 G C \ r. Thus is invertible. Consider Since Y^""^ 

and have the same jumps, this function has no jump on P. Therefore 

its entries are entire functions. Moreover, from the asymptotic condition at 



7r„+i(2:) = {z- bn)TTn - a„7r„_i(2:) 



(2.5) 




n 
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infinity it follows by Liouville's theorem that there exist constants a„, 6„ and 
c„ such that 

y (n+l) (^) (y W (^)) -1 = /"^ - -flnK-A _ (2.6) 



So 



\ Cr), / 



and this is in fact the difference equation in the Lax pair in ^H]- The 11-entry 
of equation (|2I7|) leads to (Q. □ 

The recurrence coefficients a„ and 6„ can be expressed directly in terms 
of the solution Y^^^ of the RH problem. This follows as in ^ Section 3.2]. 
For convenience of the reader, and since the precise form ()2.1U|) for the recur- 
rence coefficient 6„ is not given in [9], we include the proof of the following 
proposition as well. 



Proposition 2.3. // we write 

Y^-\z) = (/ + r/"V^ + yt^/^' + 0{z'')) (2.8) 



as z ^ oo, then 



and 



-'2 



Yl"') . (2.10) 



)u 

Proof. As in the proof of Proposition 12. 21 we have that (Y^"^^)^^ is a polyno- 
mial of degree n — 1 with leading coefficient Hn-i- Thus from ()2.8j) we obtain 

fr/")) . (2.11) 



21 

Multiplying both sides of the identity (j2.7|) with Q^^d using (j2.8j) . we 

find 



ii + o(i/.)) j = f J" J ( / + ^ + ^ + o(. 
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as 2; oo. Taking the 12-entries on both sides, we obtain 



^ '12 /€„_! V V / 12 V /12 l^n-l ^ ^22 ' 



Thus 



and 



y[n)\ ^ ^2.12) 

12 



F^^")) -fenfr/"^) -^fr/"M =0. (2.13) 

/I2 V /12 Kn~\ V /22 

Then pTT]) and (j2T2|l lead to (jT9|l . while we obtain (j2ln|l from solving 
(jrnUl for 6„ and using pT^ . □ 

Proposition 12.31 shows that we can compute the recurrence coefficients 
from the RH problem for F*^") alone. In our final proposition of this section 
we show that if we compute a„ and 6„ as in ()2.9p and ()2.10|) and if a„ 7^ 0, 
then these are indeed the recurrence coefficients. 

Proposition 2.4. Suppose that the RH problem ()2.3|) has a solution F^"^ 
with expansion (j2.8|) . Let an and bn be given by i]2.9\\ and ()2.10j) . //a„ 7^ 0, 
then the monic orthogonal polynomials Hn+i, T^n and TCn-i exist and they are 
connected by the three-term recurrence (|2.5|) . 

Proof. Since an 7^ 0, we get from ()2.9|) that 

(f/"))^^^0. (2.14) 

This means that the 21-entry of y is a polynomial of exact degree n — 1 and 
therefore the orthogonal polynomial tt^-i exists. 
Define 

On = Kn-l/an = f(r/"))l2)"' (2.15) 



and use (EHj) to define y("+i). Then is analytic in C \ T and satisfies 

the jump relation of the RH problem ()2.3|) . We show that it satisfies the 
asymptotic condition with n replaced by ri + 1. We have 

as 2 — > oo. Collecting terms according to powers of z, we find 
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there is one term with which comes with coefficient ^ J = 
O. 

the terms with z come with coefficients 

Oyi 1 vo V V cn vo V vo 

and this is O because of (QUI) and (imi . 

the constant terms are 

1 0\ A 0\ -Ctn/ft^n-A ^(n) /O 0\ A 0\ („) /O 



oyio oy'^ic^ lo ij + vo oj^^ lo i 



oy^lo c.(y/"))i2 / Vo 



;i2 



and using (Q, (EHHI), (EH), and (ITTHll . we find that this is the 
identity matrix I. 

It follows from all this that 

'^n+l Q 



= {l + 0{z'^)) 



and so solves the RH problem (j2.3|) with n replaced by n + Then the 

monic orthogonal polynomial iTn+i exists and the recurrence (j2.5j) holds. □ 

In our main results we are interested in the asymptotic behavior of a„^„(t) 
and bn,nit). As a result of Propositions 12.31 and 12.41 it suffices to take N = n 
in (j2.3|) and do a steepest descent analysis for this RH problem. 



3 Steepest descent analysis in the regular case 
and the proof of Theorem 11.11 

In the regular case — l/12<t<0 with t fixed, the steepest descent analysis 
for the RH problem ()2.3|) with N = n can be done in a way very similar to 
the classical steepest analysis for t > as in [T^ lllj. 
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In the first transformation we use tlie continuation of the equilibrium mea- 
sure that we discussed in Section 1 for t > 0, see (|3.ip . The local parametrix 
at the end-points is constructed out of Airy functions, as in the case t > 0. 
Therefore the analysis will not be carried out in full detail. We only point 
out how one gets into a standard situation, by deforming the contour and 
defining the correct equilibrium measure. Then we only sketch the rest of 
the analysis and state the result for the recurrence coefficients. 



The equilibrium measure 

The first transformation is based on the equilibrium measure fit defined in 
(II. So we have 

^ = ^(^'-c^?)v/^F^, forx6[-Q,Q], (3.1) 

with Ct and dt given by ()1.5|1 and ()1.7|) . Since —1/12 < t < 0, we have that 
dt is real and < q < dt- Also note that the density ()3.1|) is positive for 
X G (ct, Ct). The measure fit minimizes Ivtif^), ^^e ()1.4|1 . among all probability 
measures supported on [— 

Let gt : C \ (— oo, q] — C be the 5^- function defined by 

/Ct 
\og{z - x) d/it(x), (3.2) 
■Ct 

where for each x G (— Q, q) we choose the principal branch of the logarithm 
log(2; — x). Then gt can be represented in the following way 

9ti^) = A [\'^' - dt)is' - c\fl^ d. + \Vt{z) + /*/2, (3.3) 



for all z G C\ (— cx), q] and for some constant It- Define 0^ : C\ (— 00, q] C 
by 

0,(.) = -^£(.2-d?)(.^-c?)V^d.. (3.4) 

In the first transformations the jumps for the transformed RH problem are 
expressed in terms of 

Note that the integrand in ()3.4j) is negative for s G [ct^dt). Therefore 

0t(z)<O, for z G (Q,rfi). (3.5) 
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By symmetry we have 

(j)t+{z) + (j)t-{z) < 0, for z e {-dt, -ct). (3.6) 
Then it follows from (jSSl), dSSl) and (jSSD that 



gt+{z) +gt^{z) - Vt{z) = lt, z e (-q,q), 
9t+{z) + gt-{z) - Vt{z) <lt, z e {-dt,-Ct) U {ct,dt). 



The transformation Y T 



(3.7) 



The transformation Y T normalizes the condition at infinity and serves 
as a first step to get the jumps close to the identity. Define T : C \ F by 

for all 2; G C \ r. Then T satisfies the following RH problem 
T{z) is analytic in C \ F 

/Q-n{gt + (z)-9t-{z)) Q,U\Q-n(Vt(z)+lt-gt+iz)-9t-{z))\ 

T+iz)=T_iz)\^ e"(f*+(^)-^*-(^)) J' 

T{z) = I + 0{l/z), z^oo. 

(3.9) 

Since gt is analytic in C \ (— c>o, q], one might expect a jump on (—00, —dt). 
However, since gt-^-{z) — gt-{z) = 27ri for z G {— 00, —dt), the function 
Qn{gt_f_{z)-gt_iz)) -g ^^alytic for z G {—00, —dt) and therefore there is no jump 
for S on (— cxD, —dt). 

Now from ()3.3p and ()3.4p it follows that the jump matrix can be rewritten 

as 

/ /1 „J ~\„2nd>t(z^\ 

z G r \ 

T^{z)=T.{z){ ^ y ^ (3.10) 

z G (-Ct,Ct). 

For z G (— Q, q) we get from (|3.4p that (t)t±{z) G iM. The diagonal entries of 
the jump matrix on (— Q, q) are therefore oscillating functions. By (|3.5j) and 
fl3.6|) the jump matrix converges pointwise to the identity for z G {—dt, —Ct) U 
{ct,dt). We first have to specify the precise deformation of the contour to 
control the behavior of the jump matrices on Fi, . . . , F4. We will deform the 
contour such that the jump matrix converges pointwise to the identity for 
z G F \ [-Ct,Ct]. 
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Figure 4: The shaded region is the region where Re 0j > in the regular 
case. The curves are steepest descent curves for (pf The figure corresponds 
to the value t = -1/24. 

Deformation of the contour 

In view of the jump matrix (j3.1(Jj) on [— Q, q] we want that the real part of (pt 
is negative on r\ [— q, q]. The shaded region in Figure |^ is the region where 
Re 4>t{z) > 0. In the white unshaded region we have Re (pt{z) < 0. So we 
want that each Tj stays in the white region. One possible way to define F^- is 
to take the steepest descent curve for 0j through ±dt in the k-th quadrant, 
where we have that Im (pt is constant, and Re (j)t{z) — > — oo as 2; — oo along 
the steepest descent curve, see Figure El 

Summary of the rest of the steepest descent analysis 

In the remaining part of the steepest descent analysis we open a lens around 
[— Q, Ct] to turn the oscillating jumps on [— q, q] into a constant jump at 
[— Q, Ct] and a jump that converges pointwise to the identity on the upper 
and lower lips of the lens. In Figure the contour is shown after opening of 
the lens. This lens can be defined in a completely similar way as in the case 
t > 0. After this, we construct a local solution, called parametrix, around 
the endpoints iQ. But locally around ±q, our situation does not differ from 
the case t > 0: the equilibrium measure vanishes with an exponent 1/2 at 
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Figure 5: Opening of the lenses in the regular case. The lens around [— Q, q] 
is contained in the shaded region of Figure |3] 

±Q and the contours together with the jump matrices locally have exactly 
the same shape as in the case t > 0. So the parametrix is constructed by 
means of Airy functions as in [IT] . 

Proof of Theorem 11.11 

The result of the steepest descent analysis, which we will not perform in 
more detail here, is that the RH problem for Y^"^^ is solvable if n is large 
enough, so that the coefficients an,n{t) and bn,n(t) can be computed in terms 
of Y^"^^ as in Proposition 1231 It also follows that the asympotic formula ()1.8p 
for the recurrence coefficient a„^„(t) which holds for t > continues to hold 
for —1/12 < t < 0. Then a„^„(t) 7^ for n large enough, so that a„^„(t) 
and hn,n{t) are indeed the recurrence coefficients by Proposition 12.41 As in 
[TTj we also find a full asymptotic expansion for an,nii) in powers of n~^. In 
addition, the recurrence coefficient bn,nif) is exponentially small. 
This completes the proof of Theorem 11.11 
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4 Steepest descent analysis in the critical case 
and the proof of Theorem 11.21 

The rest of the paper is devoted to the critical case. We take = n in 
(12 .31) and analyze the RH problem in the double scaling limit. First we let 
t be a fixed value close to —1/12, which could be less than —1/12. Later 
we let if: —1/12 and n oo simultaneously and show that the Painleve I 
equation appears. 

In the critical case the measure ^cr is given by 

du (8 - x^)^/^ 

dx ^ ^ 247r ' ^ ^ 

for X G [-V8, V8]. The measure vanishes with an exponent 3/2 at x = ±v8- 
A consequence of this phenomenon is that the parametrix with Airy functions 
which we could use in the regular case fails to work in this case. 



4.1 Deformation of the contour 

As in the regular case we begin by deforming the contour F. The deformed 
contour is constructed in terms of (per and will therefore not depend on t. 
Here (per is 

M^) = 0-i/i2(^) = j^r is' - Sf/'ds (4.2) 

J Ccr 

and Ccr = V8- 

The shaded region in Figure IHl represents all 2; G C for which Re (pcr{z) > 
0. The white region represents all 2; G C for which Re (pcr{z) < 0. As in the 
regular case, we want the F^ to be contained in the white region. Comparing 
with Figure m we see that the shaded parts meet at ±a/8, which is a direct 
consequence of the exponent 3/2 vanishing of the measure at ±Ccr = ±V8. 
Hence a major difference with the regular case is the fact that the contours 
have to leave the real axis at the point ±Ccr = ±-\/8 to continue in the 
complex plane, whereas in the regular case we can deform the contour such 
that it contains the interval [—dt, dt] with dt > Ct. 

Now define Fj, j = 1, . . . , 4 as the steepest descent curve for (per through 
±V^ in the j-th quadrant, where we have that Im (per is constant, and 
Re (per{z) —00 as z ^ 00 along the steepest descent curves. We also 
put Fq = [—\/8, -\/8], see Figure El 
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Figure 6: The shaded region is the region where Re 0cr > 0. The curves 
Fi, . . . , r4 are the steepest descent curves for 0cr that contain ±Ccr = ±V^. 

4.2 Modified equilibrium problem 

An important feature of our analysis is the use of a modified equilibrium 
problem. This is inspired on the papers OE], see also where the authors 
used a modification of a similar equilibrium problem to handle a double 
scaling limit. 

Recall the energy functional Ivt as defined in (jl.4j) . Instead of minimiz- 
ing among Borel probability measures on M, we consider the following 
extremal problem. Minimize 

IvA'^) = jj log d'^(^) d'^(^) + / ^*(^) d'^(^) (4-3) 

among all signed measures v such that 

jdu = l, and supp (z/) C [-^8, Vs]. (4.4) 

From general arguments of potential theory ^2^, there exists a unique 
minimizer Ut of Ivf We will use this measure in the steepest descent analysis. 
Note that its properties differ on two crucial points from the equilibrium 
measure fit as defined in p.ip . It is forced to have support within [— v^, a/8], 
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and it is not necessarily positive everywhere. Only ior t — tcr — —1/12, we 

have Ut = fif. 

The use of this modified equilibrium measure is not without a price. If 
t > —1/12 the density of the equilibrium measure is negative near the end- 
points ±-\/8, which causes a problem when we open the lens. Instead of 
exponential decay, we get exponential growth in a neighborhood of ±-\/8. 
This neighborhood becomes smaller when t tends to —1/12 and will eventu- 
ally fall inside the region where we are going to make a special parametrix 
anyway. On the other hand, if t < —1/12 the measure ut is positive and the 
opening of the lens causes no problem, but we still get exponential growth in 
a neighborhood of the endpoints ±-\/8 which now follows from the behavior 
of its associated function on the contours Tj, j = 1, . . . , 4. However, when t 
tends to —1/12, this neighborhood will eventually fall inside the region again 
where we make the special parametrix. 

We will derive an explicit expression for the minimizer. The Euler- 
Lagrange equations for the minimization problem yield that for some con- 
stant If we have 



2 / log -^dutiy) + Vt{x) +lt = 

J \x-y\ 

for X E [— \/8, VS]- Taking a derivative with respect to x, we get a singular 
integral equation for the density Vf of i/f, 

2PV ^^dy = V^{x) = V;,{x) + {t + l/12)x^ for x e {-Vs, Vs). 



Thus 



where 



^(X) = Vt{x) = V,r{x) + {t+ 1/12)V°{X) 

dx 



and v°{x) is such that 



and 



v°{x)dx - (4.5) 

-^dy = xy2, X e (-VS, VS). (4.6) 
-VEX-y 
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We determine v° from (j4.5|) and (j4.6|) . 
Lemma 4.1. VFe have that v° is given by 



Proof. Define 



8 -I- At^ 

v%x) = 7^4 , X E i-Vs, V8). (4.7) 



S _|_ 4^2 _ 4 



where the branch of the square root is chosen which is positive for real z > 
Vs. Then, 

^^^^f = -ih-{y) - K{y)), y e (-v^, Vs), 

which imphes that 

8 + 4y ^ 1^ f ^ 

V8 27rv/8^ 27ri 



and 



PvT 't^-"' d, = ^/^d.. .e( 

J-^ 27rv/8-y2(a; _ 27ri Jc x - z 

where C is a closed contour around [— -\/8, Vs] with counterclockwise orien- 
tation. Since 

S _|_ 4^2 _ 4 

which we can check by straightforward calculation, we find by the residue 
theorem 



/ 1^ — /n = - Res h{z) = 0, 



y_v^ 27rA/8 - y^{x - y) ^=<« a; - z 

Thus v° defined by ()4.7|) satisfies ()4.5|) and ()4.6|) . This proves the lemma. □ 
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4.3 Modified gt and (pt functions 

We use gt to denote the ^f-function associated with vt- So : C \ (— oo, \/8] 
is defined by 

gt{z) = J \og{z - x)dut{x) (4.8) 

and it satisfies 

gt+{z) + gt-{z) - Vt{z) = k, ze {-Vs, Vs) (4.9) 

and 

gt{z) = log z + 0{l/z), z^oo. (4.10) 
We also define 0j : C \ (— oo, ^/S] by 

Mz) = Mz) + it + 1/I2)(j)°{z) (4.11) 

where 



Then 



M^) = -IVtiz) + gt{z) - ht, zeC\ (-00, VS], (4.13) 

which will be used in the first transformation. 

Observe that gt and (pt are different from the functions with the same 
name that were used in Section El in the steepest descent analysis in the 
regular case. Since we will not use the results of Section 01 in what follows, 
we trust that this will not cause any confusion. 

4.4 The transformations Y T S 

Define T : C \ T by 

T{z) =e-"'*"^/V(2)e-"^*(")"«/2e'''*"^/^ (4.14) 
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Figure 7: The contour S5 that arises after opening of the lens around 
[— Ccr, Ccr] where = V^- 

for all z G C \ r. Then T satisfies the following RH problem 
T{z) is analytic in C \ F 

T+iz) = T^z) "^^'^^ j, zGF,, J = 1,2,3,4, 

[ T{z) = I + 0{l/z), z^oo. 

(4.15) 

Since (l)t±{z) G iK for z G (— V^, v^), the jump matrix is oscillatory on 
this interval. In the second transformation we open a lens. This lens will not 
depend on t. Add two curves to the contour such that both contour go from 
-^8 to ^8, one below and the other one above this interval, such that 

Re (pcriz) > 0, (4.16) 

holds on the lips of the lens. Figure IHl indicates that this can be done. It 
leads to a contour as in Figure d where the upper and lower lips of the 
lens are in the shaded region of Figure IHl Define 
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S{z) = T{z) ^ ^-2n<j)t{z) ^ ' upper part of the lense, 

S{z) = T{z) (^^-2n<t>t{z) ^ 5 the lower part of the lense, ^^-^'^^ 

S{z) = T{z), elsewhere. 

Then S satisfies the following RH problem on the contour consisting of 
r and the upper and lower lips of the lense. 

S{z) is analytic in C \ S5 

1 a(2)e2"'^'(^) 



S+{z) = S^{z) ( oL^(^) *^ ) , on the lips of the lense, (4.18) 



S+iz) = S^iz) j, zET„ j = l,2,3,4, 

1 0\ 

2n((>t{z) ]^ I ) 

^ S{z) =I + 0{l/z), z^oo. 

The contour is constructed such that Re 4>cr{z) < on Tj, j = 1,2,3,4 
and Re 4>cr{.z) > on the upper and lower lips of the lens. So the jump 
matrices in the RH problem for S would have the correct decay properties 
in case t = —1/12, since then 0^ = (per- However, if t 7^ —1/12, then (pt is 
different from (per, and Re (pt{z) has an incorrect sign on some of the contours. 
In Figure |H1 the contour is shown near a/8 . The shaded regions are the 
regions where Re 0t > 0. The left picture shows the behavior of Re (pt for 
t < —1/12 and the right picture the behavior for t > — 1/12. In the left 
picture we see that a/8 is inside the shaded region, whereas it is inside the 
white region in the right picture. Therefore, Re (pt has the wrong sign on 
(small) parts of Tj, j = 1, 2, 3, 4 in the left picture. In the right picture it has 
the wrong sign on (small) parts of the lips of the lens. However, the following 
proposition states that this behavior only occurs in small neighborhoods of 
±78. 

Proposition 4.2. Let U and U be neighborhoods of \/8 and —\/8, respec- 
tively. Then there exist 6 > and e > such that for all t E with 
\t + 1/12| < 6, we have Ke (pt > e on the upper and lower lips of the lens 
outside UUU and Re (pt{z) < -^l^l"^ on Tj \{UUU) for j = 1,2,3,4. 

Proof. There exists an e > such that Re (per > e on the lips of the lens 
outside U and Re (pcr{z) < —e\z\^ on the contours Tj outside U. The latter 
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Figure 8: Behavior of Re (pt in a neighborhood of Ccr = V8 for t < —1/12 
(left) and t > —1/12 (right). The shaded region is where Re 0t > 0. The 
real part of (pt has the wrong sign on some parts of S5. 

inequality holds since (pcri^) = (1/96)^"^ + 0{z^) as z — 00 by (j4.2j) and Tj 
is a steepest descent curve for (per- Since 

4>t = (per +{t + 1/12)0°. 

and (p°{z) = —{l/8)z^ + 0{z^) as 2 ^ oo, see ()4.12|) . it follows by continuity 
that the same inequalities hold for (pt if t is sufficiently close to —1/12. □ 
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4.5 The parametrix M away from the endpoints 

The 2 X 2 matrix valued function M : C \ [-v^, v^] ^ C^^^ defined by 

is a solution of the RH problem 

{M is analytic in C \ [-v^, ^8], 
M+{z) = M.{z){j^^ J^, ^G(-v^,y8), (4.20) 

M{z) = I + 0{l/z), z^oo. 

Thus M has the same jumps as S on the interval (— -\/8, Vs). We will use 
M as a parametrix for S away from the endpoints. That is, we use M as an 
approximation for S away from ±a/8. The approximation will be good if n 
is large, and t is sufficiently close to —1/12. The approximation will not be 
good near the endpoints. There we construct a local parametrix with the aid 
of the RH problem for Painleve I which was already given in Section ^ but 
which will be discussed in more detail next. 

4.6 The RH problem for Painleve I 

In the literature one can find several equivalent versions of the RH problem 
for Painlevee I [201 1211 ■ Here we start from the RH problem as formulated 
by Kapaev [22] . 

Let Exj, be the collection of rays arg(C) = 2A;7r/5 for k = —2, —1,0, 1,2 
together with the negative real axis arg(C) = tt with orientation as in Figure 
[3 Let v4o(C) for C G be as indicated in Figure [HI that is, 

arg(C) = 2A;7r/5, A; = -1, 1 

arg(C) = 4A;7r/5, A; = -2, 0,2 (4.21) 

arg(C) = vr, 

where the Stokes multipliers Sk are complex numbers such that 

1 + SkSk+i = -iSfc+3, = Sk, ke Z. (4.22) 



^(0 = <^ 
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Figure 9: The contour Sij, and the jump matrices Aq for the RH problem 
for 

Consider the RH problem for a 2 x 2 matrix valued function "^^^^ depend- 
ing on a parameter x, 

x) is analytic in C \ S^, 
M/f (C;a;) = ^^'\c,x)Ao{C), ( e E^, ^^ ^S) 

where 6 is given as before by (IT^ . Then has an expansion of the form 

(4.24) 

where TC and ?/ depend on x. The function y = y{x) is a solution of the 
Painleve I equation ()1.12|) and TC{x) = — 2y^{x) — yx is the Hamil- 

tonian. Note that TC'{x) = —y{x). To every set of Stokes multiphers Sk sat- 
isfying (j4.22j) . there corresponds a unique solution of (jl.l2|) and vice versa. 
The RH problem for '^^^^ has a solution if and only if x is not a pole of y. 

For our purposes we need a special choice of Stokes multipliers. We take 
So = and si = ia. This determines the other Stokes multipliers by ()4.22p 
and it follows that 

So = 0, si = ia, s_i = i(l — a), S2 = s_2 = i. (4-25) 
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We will denote this special solution by x, a). The corresponding so- 

lution of the Painleve I equation and its Hamiltonian will be denoted by ya 
and Ha- 

Before we proceed with the parametrix we first modify the RH problem 
for \E'*^*^) by defining 

vi/(C;x,a) = vI/W(C;x,a)(^J (4.26) 

and we reverse the orientation on the contours arg^ = vr and arg^ = ±47r/5. 
Then the RH problem for reads 

\E'(-; X, a) is analytic in C \ 

^^+(C; X, a) = ^4C, X, a)A{C), C e E^,, 

vI/(C;x,a) = ^(j 7)(/ + ^ + ^ + 0(^))e^(^-)-, C - oo. 

(4.27) 

where the jumps A{Q are as indicated in Figure El in Section [T] and 



where we use 



1\ fO -i\ fl 



to denote the Pauli spin matrices. 

We note that (^~'^3/4v|>(o)g-6'o-3 j^^^g ^^jj asymptotic series in powers of 

C"^''^, see |2m I28j. Therefore \E' has also an asymptotic series 

v^(C; «) ~ ^ (1 ~^ (^^ + E «)^"'^'] (4.30) 

for C — > oo. For our purposes it is sufficient to work with the terms up to 
order 0{(~^). It is possible to derive a full asymptotic expansion for the 
recurrence coefficients if one works with the full expansion ()4.30p . 
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4.7 Construction of the local parametrix 

We will define local parametrices P and P for the RH problem for S around 
the endpoints ±-\/8- The jumps for S are given in ()4.18|) and we want that 
P has the same jumps in a neighborhood of ^/S. In addition we want that 
P matches with M that we constructed in ()4.19j) . So in a disk U around 
we want that P satisfies the following RH problem. 

P{z) is analytic in f/ \ S5 

on the lips of the lens inside U, 

p+iz) = p^{z) i^^^ qJ, zEi-Vs,V8)nu, 

P{z) = M{z)(I + 0{n~^^^)), as n — > 00, uniformly for z E dU. 

(4.31) 

The construction of P is based on the following observation. 
Lemma 4.3. Let P -.UXJ^s ^ C^^^^ be defined by 

P{z) = E{z)^{n^/^f{zy, n^/^utiz), a)e-"'^'(^)"3 (4.32) 

where 

(1) '^{■■,x,a) is the solution of the RH problem (j4.27|) that is associated 
with 

(2) f is a conformal map from U to a neighborhood of such that S5 fl t/ 
is mapped to (part of) the contour S^; 

(3) Ut : U ^ C is analytic and n^^^Ut{U) does not contain any poles of ya,' 

(4) E -.U ^ C2><2 IS analytic. 

Then P is analytic m f/ \ S5 and satisfies the jump conditions in fl4.31|) . 



33 



Proof. If we put P = Pe^'t'*^^ then the jump properties for P in (|4.3ip trans- 
late into 



P+{z)=P^{z) 
P4z) = P4z) 

h{z) = p-{z) 

Pjz) = P(z) 



1 a 

1 

1 1 - a 

1 

1 0\ 

1 1 ' 

l\ 

-1 ' 



zeTinu, 

z eTiHU, 

on the lips of the lens inside U, 

z e {-Vs, Vs) n u. 



(4.33) 



Here we used the fact that (pt is analytic in C \ {—oo, ^/S] and (pt^ + (f)t_ = 
on {-y/8,y/8). The jump matrices in ()4.33|1 are exactly the same as the ones 
that appear in the RH problem for see Figure El except that these jumps 
are on different contours. Since / provides a conformal map between the 
respective contours, it is then immediate that '^{n'^^^f{z);x,a)e~'^'^*^^^'^''^ has 
the required jump properties on E5 fl U, for every x that is not a pole of i/a- 
Since the jumps do not depend on x, the jump properties will not be affected 
if we let X = n'^/^Ut{z) where Ut is analytic and n^^^Ut{z) is not a pole of 
ya for every z & U. Finally, we note that the multiplication on the left by 
an analytic factor E{z) does not change the jumps either, so that ()4.32|) has 
indeed the required jump properties on fl f/. □ 

What remains to be shown is that we can find /, Ut and E with the 
properties stated in Lemma f4.3l such that the matching condition 



P{z) = M{z){I + 0{n-^l^)), as n ^ cx), uniformly for z G dU, (4.34) 



is satisfied as well. In the next lemma we will show that this is indeed satisfied 
if we define 

'5 , , 



2/5 



(4.35) 



(4.36) 



and 

E{z) = M{z) 



{n^">f{z)Y^/^ (I -i 
V2 V i 



(4.37) 
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Recall that (per is given by (j4.2j) . 0t is given by (j4.11|) . and M is given by 
(|4.19p . The branches of the fractional exponents in (j4.35|) - (|4.37|) are all taken 
to be positive for z — y/S real and positive and sufficiently small. 

Lemma 4.4. Assume that x is not a pole of and let 

t = tcr- ClXn-^/^ ci = 2-9/^3-^/^ (4.38) 

Let f , Ut, and E he as in (j4.35j) - (j4.37j) . Then there is a disk U around -\/8, 
depending only on x and a hut not on n, such that the properties (2)-(4) 
stated in Lemma 14.31 are satisfied and such the matching condition ()4.34p 
holds. 

Proof. We begin with the properties of /. It follows from this definition ()4.2p 
of 0cr that 

27/4 

(Pcriz) = —{Z - VSf/^l + Kr{z)) (4.39) 

with an analytic function in a neighborhood of -\/8 with /icr("\/8) = 0. 
Thus (|4.35|) defines a conformal map ( = f{z) from a small enough disk Ui 
around a/8 to a neighborhood of C = 0. Then f{z) is real for real z G f/i and 
since by fl4.39|) and ()4.35|) we have 

f{z) = 2-i/i03-2/5(z - 78) + 0{iz - Vsf), (4.40) 

as 2; — >• v^, we see that / maps {—Vs, Vs) (1 Ui to a part of the negative 
real axis. 

Since Fi and r4 are defined as the steepest descent curves of (per that 
start from v^, we have that (per is real and negative on Fi and Then it 
follows from (j4.35j) that / maps these contours onto the rays a.Tg( = 27r/5 
and argC = — 27r/5, respectively. 

Finally, to have all the mapping properties stated in (2) in Lemma 14. 3| 
we want that / maps the upper and lower lips of the lens that are inside Ui 
to parts of the rays a.Tg( = ±47r/5. Here we use the additional freedom we 
have in choosing the exact location of the lips of the lense. So far we have 
only specified that they should be in the region where Re (per > 0. Now we 
require in addition that in a neighborhood of we want the lips to be such 
that a,Tg(per{z) = ±2n for z on the upper and lower lips near v^. We can 
clearly impose this extra requirement. Note that we take a.Tg(per{z) so that 
it is continuous for 2; G C \ {—00, ^/8] and has the value for 2; — a/8 real and 
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positive. Then by (j4.35p we have that arg/(z) = ±47r/5 for z & Ui on the 
hps of the lense, and thus all the properties stated in (2) of Lemma 14.31 are 
satisfied. 

Now we turn to the properties of Ut- Because of fl4.1ip and ()4.36p we can 
write 

utiz) = {t + l/12)u°{z) (4.41) 

where 

From KT^ it follows that 

(j)°{z) = -2^/^3 {z - Vsy/^il + h°{z)), (4.43) 

where h° is analytic in a neig hborhood of ^8 with h°{V8) = 0. Then by 
fl4.39p . ()4.42|) . and ()4.43|) we find that u° is analytic in a neighborhood of \/8 
and 

u°{z) = -29/^3^/^ (^1 + Oiz - VS)) (4.44) 

for z ^ VS, so that u°{V8) = -c^\ Combining (OH|l . (OT|) . and (I01|l . 
we see that n^^^Ut{V8) = x. Since x is not a pole of ?/q, we can then find a 
disk U2 around a/8 with f/2 C f/i, such that n^^^ut{z) is not a pole of |/q, for 
every z G ?72- Since n^^^Ut = —Cixu° depends only on x, the disk U2 depends 
only on x. So the properties stated in (3) of Lemma [4.31 are satisfied. 

Since M is given by ()4.19p we find from ()4.37|) that 

n^/'f{z){z+ 




which implies that E is analytic in a disk U C U2 around a/8, as f{z) has a 
simple pole at z = a/8, see ()4.4()j) . This shows that property (4) of Lemma 
14.31 is satisfied. 

We finally show that the matching condition ()4.34|) is satisfied. Since \l/ 
has asymptotics ()1.2H1 . we find from ()4.32|1 that uniformly for z e dU, 

X exp {{e{n^/'f{z), n'/'ut{z)) - n<f)t{z))a^) (4.45) 
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as n ^ oo. Since 9{C:^) = IC^''^ + xC^^'^, we see that 



ncptiz) (4.46) 



where in the last equality we used the formulas (j4.35p and (j4.36|) for / and 
Ut- Hence we can forget about the exponential factor in (j4.45p . Because of 
the definition ()4.37p of E we see that ()4.45p leads to the matching condition 

This completes the proof of Lemma 14 .41 □ 

Remark 4.5. The neighborhood U in Lemma 14.41 depends on x, but an 
inspection of the proof shows that that we can take U independent of x if 
X is allowed to vary in a compact set that does not contain any poles of ya- 
Then also the matching condition ()4.34|1 is uniformly valid for such x. 

Parametrix around — 

The local parametrix P in a disk U around the other endpoint — should 
satisfy the following RH problem. 

P{z) is analytic in [/ \ S5 

/l fl - /?)e2"'^*(^)\ 
P^{z) = P4z) 1^^ j, zeT.nu, 

p+(^) = p_(^)(^J ^\ y ZGLgnf/, 

P+(2;) = P-{z) ^^_2n0t(2) , on the lips of the lens inside U, 

p+(^) = p_(^)(^_°^ J^, zei-V8,V8)nu, 

, P{z) = M{z){I + 0(n"^/^)), as n — >• 00, uniformly for z E dU. 

(4.47) 

It turns out that P can be expressed directly in terms of P. Let us write 
P = Pq, to emphasize that the solution of ()4.3H) depends on a. Then it 
follows that 

P{z) = a^Pp{-z)a; (4.48) 

solves ()4.47|) on the neighborhood U = —Up of — \/8. Indeed, the identity 
()4.48|) is an immediate consequence of the RH problems ()4.3H) and ()4.47|) . 
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and the facts that M{—z) = a3M{z)a3 and (f)t{—z) = (pti^) ± vri. It is well- 
defined if X is not a pole of y^. 

Since f,Ut, and E do not depend on a, it follows from ()4.32j) and ()4.48p 
that 

P{z) = a3E{-z)<il{n^/^f{-z); ^)e-"'^*(-")"V3 (4.49) 

for z G U. 

Matching condition 

We obtain a more precise matching condition than ()4.34|) if we use the asymp- 
totic series (j4.3Up in (j4.32j) . Then it follows as in the proof of Lemma 14.41 
that 

P{z) ~ M{z) (^I + fl ^^kin'^'utiz), «)/(^)-'=/2^-'=/5 j (4.50) 
as n — » oo, uniformly for z G dU, and 

P{z) ~ Miz)a, (^^ + £ ^k{n'/'ut{-z),(3)f{-z)-'/'n-'/'^ a, (4.51) 

as n — ^ CX3, uniformly for z G dU, where in the last identity we also used that 
a3M{-z)a; = M{z). 

For future reference we recall that 

n'^/^utiz) = -cixu°{z) (4.52) 

where u° is an analytic function with 

u°{V8) = -q^ (4.53) 

4.8 Third transformation S R 

Now let X be fixed and assume that x is not a pole of and yj3. As before 
we let t vary with n such that 

t = -1/12 - ciCT-^/^ (4.54) 
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Figure 10: Contour Tiji for the RH problem for R 



Then the local parametrices P and P are defined in U and U and they 
satisfy ()4.31|) and ()4.47|) . We also have M as a parametrix for S away from 
the endpoints. 
Define R as 

( S{z)M-^{z), ^ G C\ (^Fut/USc.) , 
Ri^) = I Siz)p-\z), zeU\J:s, (4.55) 

[ Siz)p-\z), zeu\i:s- 

Then R satisfies a RH problem on the contour S/j that is shown in Figure 
ITUl consists of the curves S/jj for j = 1,2,3,4 which are the parts of 
the curves Tj outside the disks U and U, of the curves S/j^s and S/j g which 
are the parts of the upper and lower lips of the lens that are outside the two 
disks, and of the circles dU and dU. The orientation of Hji is as indicated in 
Figure ITUl In particular we choose clockwise orientation on the two circles. 

Lemma 4.6. With S/j as described above we have that R satisfies the RH 
problem 

R is analytic in C \ S/j 

R+iz)=R.iz)ARiz), zeT.R, i = 1,2,3,4, (4.56) 
R{z) = I + 0{l/z), z^oo, 
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where 

Ar{z) 
ARiz) 



M{z) 



1 a(2)e2"<^*(^) 
1 



M~Hz) 



z&^R,j, j = 1,2,3,4, 

(4.57) 



M{z) 



M~Hz) 



1 0^ 

g -2n<f>t{z) 

P{z)M{z)-\ z e dU, 

P{z)M{z)-\ zedU. 



Z G ^R,5 U ^R,6, 



(4.58) 

(4.59) 
(4.60) 



Proof. Note that 5* and M have the same jumps on (— v^, V^), S and P have 
the same jumps on E5 fl U, and 5" and P have the same jumps on fl [/). 
From this it follows that -R+ = -R_ on these contours, so that R is analytic 
there, and these contours do not appear in T,r. 

The jumps ()4.57|) - ()4.60|) are an easy consequence of the definition ()4.55p 
and the jumps satisfied by S, see fl4.18|) . □ 

The jump matrices (|4.57p and (j4.58|) are exponentially close to the identity 
matrix as n — 00. This follows from Proposition 14.21 since t — * —1/12 as 
n 00. 

Because of fOni) - (H3T|l and (0321) -(ESOl) we have that Ar{z), z e dU U 
dU, has an asymptotic series in powers of n~^/^ 



Ar{z) 



/ + J]^y(^)(z)n-^/^ zedU, 

k=l 

00 



(4.61) 



k=l 



where 

lyW(z) = W^'Xz) = M{z)^k{n^/''ut{z),a)M'\z)f{z)-''/\ 
and, see also (|4.48p . 

W^'\z) = a,W^'\-z)as. 
It follows that R admits an asymptotic series in powers of n~^^^. 



(4.62) 
(4.63) 
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Lemma 4.7. For large enough n, the RH problem for R has a solution, and 
R admits an asymptotic expansion 



R{z) ~ / + 5^ R^''^ (z)n-^/^ (4.64) 



k=l 



as n ^ OO, uniformly for z G C \ Sr. 

The expansion ()4.64|) is valid uniformly near infinity in the sense that for 
every K > 1 there is a constant Ck > such that 



A'-l 

" < CK\z\-^n-^^/^ (4.65) 



R{z) - I - J2 

k=l 



holds for every z with \z\ > 3. 

Proof. We already observed that the jump matrix Aji for R is close to the 
identity if n is large. On the parts Hrj we have by Proposition 14.21 that 

\\Ar{z) - I\\ < Ce-'\'\\ ze^Rj, j = l,...,6, (4.66) 

for some positive constants C and e. 

Then the lemma follows from ()4.66|) and the expansions ()4.61|) . cf. also 
the arguments used in the proof of |12| Theorem 7.81]. 

The estimate (|4.65p follows as in Lemma 8.3]. □ 

Corollary 4.8. The RH problem for Y is solvable for sufficiently large n. 

Proof. Since the transformations Y v-^ T ^ S ^ R are invertible we can 
recover Y from R. Since the RH problem for R has a solution for large enough 
n, it now also follows that the original RH problem for Y has a solution for 
large enough n. □ 

From Proposition l2. H it now follows that the monic orthogonal polynomial 
T^n^nit) exists for n large enough. In the rest of the proof we will calculate 
the coefficients an,n(^) and hn^n{i) by means of Proposition 12.31 Then we find 
that an,n{t) 7^ for large enough, so that by Proposition 12.41 the monic 
orthogonal polynomials of degrees n + 1 and n — 1 exist as well, and that 
ci-n,n{t) and bn,n{t) are the recurrence coefficients in the three-term recurrence 
relation. 
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4.9 Explicit expressions for and i?^^) 

The matrix coefficients R^''\z) that appear in the expansion ()4.64|) can be 
found by solving iteratively the following RH problems 

R^^') is analytic in C \ {dU U dU) 

Rf{z) = R^'\z) + Et, Ri'-^z)W^^Kz), z e dU, 

R^i\z) = R^^\z) + Eti R^--'\z)W^Kz), zedU, ^ ' ' 

^ R^''\z) = 0{l/z), z^oo, 

where R^^\z) = I. These are additive RH problems and therefore can be 
solved by the Sokhotskii-Plemelj formula. It turns out that the jump matrices 
have analytic continuations inside U and U with a pole at ±-\/8- This makes 
them easy to solve in an explicit way as in |2E]- 

In order to establish p.28j) and p.29j) we have to determine R^^^ and R^'^^ 
explicitly. First we show that W^^^ and W^'^^ are analytic in U \ {V^}- We 
recall that C3 = 2^/^^3'^/^, which was already defined in Theorem 11.21 

Lemma 4.9. The functions W^^"^ and VT^^) are analytic in U \ {Vs} with 
simple poles at \/8, and residues 

Res W^^\z) = -2''^cfl-La{x) ((T3 - iai) (4.68) 

and 

Res W^^\z) = IcsHlix)! + Jc32/„(a;)a2. (4.69) 
z=Vs 2 2 

Similarly, we have that W^-^'' and W^-'^'^ are analytic in U \ { — V^} with 
simple poles at —y/S, and 

Res W^^\z) = 2^'^cJ'''Hp{x) (as + iai) (4.70) 

and 

Res W^^\z) = -]-c;nl{x)I + ]-c^yp{x)a2. (4.71) 
\/8 2 2 



z 



Proof. Because of and we obtain from that 

W^'\z) = -'^^^^t^0^M{z)a,M-\z), (4.72) 
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and 



^2 



2 fiz) 2 fiz) 

It turns out to be convenient to rewrite M from (j4.19|) as 



M{z)a2M-\z). (4.73) 



M{z) = - 
and M"^ = M* as 
M-Hz) 



z — 



1/4 



(1 + ^2) + - 



1/4 



1 I z+Vs 

2 



1/4 



(/ + a2) + - 



1/4 



Then by straightforward calculations it follows that 

1/2 

(0-3 + icTi) + 



M(z)a3M-^(2) 




z — 



z + 



{I -02). 



1/2 



2+ V8 



and 



(era -lai, 

(4.74) 

M(;z)a2M-^(2) = aa. (4.75) 

By fl4.40j) we have that / has a simple zero at v^, and so we obtain from 
and (H7ill - (lt7H|l that W^^^ and VT^^) ^re analytic in U with 
simple poles at -\/8. Because of ()4.53|) and 



z — 



lim 



21/10^2/5 



C3 



(4.76) 



we find the residues (I4.68p and (j4.69|) . 

The statements (11701) and K7T\ about W^^'^ and W^'^'^ then follow from 
(gSHl), (frail , and (irail . 



□ 



Explicit expression for R^^^ 

The RH problem ()4.67p with k = 1 reads 

' is analytic in C \ {dU U 9f7) , 

R^l\z) + W^^'^{z), zedU, 



R^lHz) = R^l\z) + W^'^iz), z e dU, 

00. 



(4.77) 



[ RW{z) = 0{1/z), 
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Lemma 4.10. The solution to (|4.77|) is given by 
1 „ . 1 



Res W'^^\z) + 



^ Res W'^^\z), 

Z + V8 z=-a/8 



1 



Z — 2=v% 

— ^ Res + ^ .... , 

z — v8z=V8 2; + v8 2=-V8 

^ ResVr(^n-2) + — ^ Res - ^ G [/. 



;z e C \ (f/Uf/), 



Res W^^\z) - W^^\z), zeU, 



z 



Z 



+ i/8 2=-V8 



(4.78) 

Proo/. Let be defined by (g^HI). Since ly*^^-* is analytic in f/ with a 
simple pole at y^j and ly*^^^ is analytic in U with a simple pole at — v^, it 
is then clear that i?*^^^ satisfies all of the properties in ()4.77p . including the 
analyticity at ±\/8. □ 



If we use the explicit expressions 



and (pr7ri|l for the residues of W^'^'i 



and W'^^\ we find the following expansion of R^'^^z) as ^ — >■ oo. 
Corollary 4.11. We have that 



(1) 



+ 



+ 0{z- 



z ^ oo 



where 



(1) 



and 



d(1) 

K2 



Res W^^\z)+ Res W^^\z) 

z=V8 z=-VS 



Res W^^\z) - Res W^^\z) 



(4.79) 



(4.80) 



(4.81) 



Explicit expression for i?^^) 
The RH problem ()4.67j) for A; = 2 reads 
' jg analytic in C \ (dU U dU 



Rf{z) = R^'\z) + W^''\z) + R^l\z)W'^^\z), z e dU, 



,(2) 



R^IUz) 



[ R('Hz) = 0{l/z), 



R^_\z) + W^^\z) + R^}>{z)W^' 



z G dU, 
z oo. 



(4.82) 
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The jumps on dU and dU consist of two terms. Both terms have an analytic 
extension inside U and U with a simple pole at ±a/8. Therefore R^"^^ can be 
given in an explicit form in the same way as 

Lemma 4.12. The solution to (I4.82|) is given by 



1 

+ 



Res [W^^\z) + R'^^\z)W^^\z)\ 



Res 



Z + -\/8 2=--s/8 - 



W^^\z) + R^^\z)W'^^\z) 



zeC\{UUU), 
(4.83) 



w/izk for z G U (z E U) we have that R^'^\z) is given by ()4.83p minus the 
jump matrix for R^"^^ on dU , (on dU ). 

We can further evaluate the residues in ()4.83p using Lemma f4. 91 Indeed, 
from (lOHIl and we get 

Res [W^^\z) + R^^\z)W^^\z)\ 

= IcsHlix)! + ^csya{x)a2 - 2'/V/^na{x)Ri{V8){a3 - iai)(4.84) 

Since {a^ — iai)^ = 0, we find that in evaluating i?i(v^)(cr3 — icri) we can 
ignore terms in i?i(\/8) involving as — iai. Using ()4.fi8|l . ()4.7()j) . ()4.72|1 . and 



fl4.74|l in the formula ()4.78j) . we see that the only other terms in i?i(\/8) 
involve + iai. Sine (as + icri)(cr3 — icri) = 2(/ + (T2), the result is that 



R^'\V8)ia,-ia,] 



2-^/^/^'Hp{x) + 2-'"^no.{x) lim 



1/2^ 



V8 V /W 



2(J + 



2-^/^cf {np{x)+n^{x)){I + a2) 



where we used f|4.76|) to obtain the last equality. Subsituting this in ()4.84|) . 
we find 

Res \W^^\z) + R^^\z)W^^\ 
z=Vs 



C3 {nl{x)I + ya{x)a2 - Ha{x) {HM + 'Hp{x)) (/ + (T2)) 



2 

2^3 {ya{x)a2 - HaHpI - Ha{Ha{x) + H(3{x))a2) ■ 



(4.85) 
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Similarly, we have 



Res 

;=-V8 



Mx)a2 + n^HfsI - np{n^{x) + np{x))a2) . (4.86) 



If we use ()4.85|1 and ()4.86|1 in ()4.83|) and expand around z = oo, we obtain 
the following corollary. 

Corollary 4.13. We have that 

p(2) d(2) 

i?(2)(^) = :^ + ^ + 0(;,-3)^ (4.87) 
z z^ 

where 

R? = ivaix) + y^ix) - (Haix) + Hfsix))^) a2 (4.88) 

and 

= y2c3 ((i/„(x) - y^(x))(T2 - 2n^{x)rLp{x)I - (nlix) - nlix))a2) . 

(4.89) 

4.10 Recurrence coefficients in terms of R 



We start from the formulas ()2.9p and ()2.10p that express a„,„(t) and hn^nif) 
in terms of the solution Y of the RH problem ()2.3|) . Following the effects of 
the transformations Y ^ T ^ S ^ R we obtain the following expressions 
for an,n{t) and bn,n{t) in terms of the coefficients Ri and R2 in the expansion 

R(z)=I+^ + ^ + 0{z-^), z^oo. (4.90) 

z z^ 

Lemma 4.14. We have that 

anAi) = ((^1)21 - iV2) ((i?i)i2 + iv^) (4.91) 

and 

h m _ (-^1)11 ~ 2~^/^i(i?2)l2 Q^v 
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Proof. Since Ut is a measure with an even density, we have 

gt{z) = j \og{z - x) dut{x) 

= \og{z) - \ j ^ ^^t{x) - ^ y dz/j(a;) + 0{z-^) 
= \og{z) - cz'"^ + 0{z~^), z oo, 
for some constant c. Then 

SO that we get from (Hmi and (IT^ that 



e 



= / + — + ^ + 2^00, (4.93) 

z z^ 

with Ti = e-nha3/2Y^Qnka3/2 g^^^ j.^ ^ Q-nka3/2 ^ ^^^^^ gn«ia3/2_ rj.^^^ 

find after some straightforward calculations 

a„,„(t) = (^1)21(^1)12 = (Ti)2i(Ti)i2 (4.94) 

and 

bn,n{t) = 7^ - (>^l)22 = - (T,),,. (4.95) 

Since by ()4.17p and ()4.55p we have that T = S = RM outside the disks 
and the lens, and M has the expansion 



Miz) = / + ^ + ^ + Oiz--') = I-^a, + -J + Oiz-% z - 00. 

Z Z'^ Z Z'^ 



we find from (IOn|l and that 

T^ = Ri + Mi = Ri- V2<J2, T2 = i?2 + RiMi + M2 = i?2 - V2i?icx2 + /. 

(4.96) 

Inserting ^^T^ into and we arrive at (jCT]) and pH^ . □ 
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4.11 Proof of Theorem O 

Finally, we are ready for the proof of Theorem 11.21 

Proof. For each k > 1 we have that R^'^^ has a Laurent expansion at infinity 



z oo. 



Because of ()4.fi5j) and ()4.90|) we then get that Ri and R2 have asymptotic 
expansions in powers of n~^^^, 



R,--Y.^fn-''\ J = 1,2. (4.97) 
fe=i 



Using ()4.97|) in ()4.91|) and ()4.92|) we find that a„„(t) and bn,n(t) have asymp- 
totic expansions in powers of n~^^^. 
Inserting 

i?i = + i?f^n-2/5 + Oin-^/'), (4.98) 

into pin]) yields 

a„,.(t) = 2 + v^i - K^"* 
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. ((«!")., (M"),. + ((«S^'). - J) - 

(4.99) 

From the explicit expression for i?!^^ in (ICTD it follows that 

= = 2i/^cf i(7^,(x) +7^^(x)) (4.100) 

and so the order n~^^^ in (j4.99p vanishes. From the explicit expression for 
Rf^ in Km it follows that 

(i?f))^^ = - (i?f))^^ = ^C3i {y^{x) + yp{x) - (Ho^ix) + Upix))^) . 

(4.101) 
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Therefore by (j4.1UUj) and (|4.1(JH) the coefficient in the term of order n in 
flCTjl is 

= -2^/^c3{ya{x)+yf3{x)). 

This proves equation (jl.28j) in Theorem 11.21 since C2 = 2^/^03. 

The calculations for bn.n{t) are slightly more involved. Besides ()4.98p we 
also use 

i?2 = i?Wn-i/5 + + C(r2-3/5). (4.102) 

and the explicit expressions and for Rf\ j,k = 

1,2. This yields 

(i?i)n - 2-'/\R2)u = 2'/'cf{n^{x) - npix))n-'/' 

+ C3{-ya{x) + y/six) + nl{x) - nl{x))n~^"' + 0{n-^"'), 

(4.103) 



1 - 2-i/2i(i?i)i2 = 1 + 2~^l^c]l\n^{x) + np{x))n-^"' + 0(n-2/5) (4.104) 
and 

(/2i)22 = 2'"cy\n.{x) - np{x))n^'/' + 0{n-'/'). (4.105) 

Inserting (j4.10;-i|) - (j4.105|) into we obtain (fr^ . 

This completes the proof of Theorem 11.21 □ 
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